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ABSTRACT 

The ratio of ion to electron heating due to the dissipation of Alfvenic turbulence in 
astrophysical plasmas is calculated based on a cascade model for turbulence in weakly 
collisional plasmas. Conditions for validity of this model are discussed, a prescription 
for the turbulent heating is presented, and it is applied to predict turbulent heating 
in accretion disks and the interstellar medium. 
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1 INTRODUCTION 

Plasma is a ubiquitous form of matter in the universe, 
nearly always found to be magnetized and turbulent. 
In a wide variety of space and astrophysical plasmas, 
the dissipation of turbulent fluctuations makes a signif- 
icant contribution to the thermodynamic heating of the 
plasma, requiring an understanding of this behavior to 
interpret a large body of astronomical observations. Re- 
cently, there has been much vigorous activity in the study 
of turbulent heating in a wide variety of systems, includ- 
ing g alaxy clusters l|Chandranl |2004|; Dennis fe Chandranl 
2005) , the interstellar medium iMinter fc Spanglerl Il997t 




Struck fc Smith! Il999l; ISpanglerl |2003l ; IScalo fc E lmcgrccn 
20041; iLerche etal* 2007 ), accretion disks dBalbus et a" 
IQuataert 19981 ; lGruzinovlll998l; IQuataert fc Gruzinov 
. ISharma et all 120071; iRossi et all |200S|). mole cular 
clouds dPan fc Padoarj |2009l; iLemaster fc Stond 12009)) . the 
solar corona l|Vellil 120031; IVerdini et all |2009l; IChandran 



2010[ ), the solar wind llSaito et alj 120081; IValentini et al 



20081 ; [Breech et all 120091; ICranmer e t al J 120091 ; iLehe et al 
20091 ; IChandran et all l201Ch. and the magnetospher es of 
the E arth (jChaston et aT[ 20081 ) and of other planets (|Sauij 
l2004h . 

For many of these plasma environments, the dissipation 
of turbulent fluctuations occurs at length scales, in the di- 
rection parallel to the local mean magnetic field, that are 
small compared to the particle mean free paths. The result- 
ing weakl y collisional dynam ics then requires a kinetic de- 
scription ( Howes et al. 2008b J, with the dissi pation occur- 
ring via collisionless wave -particle interactions (|Howesil2008l ; 
ISchekochihin et alj 120091 ). These kinetic damping mecha- 
nisms generally transfer energy to the plasma ions and elec- 
trons at different rates. The time scale for the collisional 
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equilibration of temperature between species is often long, 
allowing the plasma ions and electrons to maintain distinct 
heating rates and temperatures. The radiation emitted from 
the hot, magnetized plasma will often depend strongly on 
the nature of the turbulent plasma heating, exerting a dom- 
inant influence on the observational signature of the object. 

Until very recently, the question of differential heating 
of the ions and electrons by plasma turbulence had not been 
widely considered, with the exception of a pioneering se- 
ries of papers b y Quataert and Gruz i nov (IQuataert! 1 19981 ; 
lGruzinovlfl99Sl ; iQuataert fc Gruzinov! 119991 ). The need for 
an accurate description of the turbulent plasma heating by 
kinetic mechanisms motivates the two primary aims of this 
Letter: (1) to calculate the relative heating of ions and elec- 
trons due to the kinetic dissipation of plasma turbulence; 
and (2) to provide a prescription for turbulent plasma heat- 
ing that may be used as a sub-grid-scale model for heating 
in simulations of astrophysical plasma turbulence. 



2 TURBULENT HEATING MODEL 

The turbulent heating presc r iption determined in this Letter 
is based on the lHowes et al.l (|2008a) model for the turbulent 
cascade of energy in a weakly collisional plasma. Whereas 
that study focused on conditions relevant to turbulence in 
the near-Earth solar wind, here we generalize to a more 
broad range of plasma parameters for general astrophysical 
applications. A brief description of the model and its under- 
ly ing assumptions follo ws — more detailed discussion is found 
in lHowes etafl (|2008al ). 

The cascade model makes three primary assumptions: 
(1) the Kolmogorov hypothesis that the energy ca scade is 
determined by local interactions (|Kolmogorovlll94ll ) ; (2) the 
turbulence maintains a sta te of critical balance at all scales 
ijGoldreich fc Sr idhar 1995); and (3) the linear kinetic damp- 
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ing rates are applicable in the nonlinea rly turbulent plasma. 
Modern theories of M HD turbulence (|Goldreich fc Sridharl 
ll995tlBoldvrevll2006h predict an anisotropic cascade of tur- 
bulent energy through wave vector space in a magnetized 
plasmaQ. The resulting sta te of critical balance is supported 
by numerical simulations dCho fe Vishniad l2000l) and ob- 
serva tions in the solar wind l|Horburv et all 1200a : IPodestal 
2009). A recent study of test particle heating in MHD turbu- 
lence supports the use of linear theory t o estimate the dam p- 
ing and heating in a turbulent plasma |Lehe et al.ll2009l ). 

Astrophysical turbulence typically supports a large in- 
ertial range at perpendicular scales larger than the ion Lar- 
mor radius k±pi < 1, leading to turbulent fluctuations 
at kj_pi > f that are highly anisotropic with fcii <C k±. 
Such anisotropic fluctuations are optimally described by 
a low- frequency expansion of kinet i c theory called gyroki- 
netics |Rutherford fc Friemanll 19681; iFrieman fc Chenlll982l ; 
iHowes et al.l 120061 : ISchekochihin et al.||2009l ). In this limit, 
all collisionless dissipation occurs via the Landau resonance, 
and the cyclotron resonance is negligible. Limitations of the 
application of gyro kinetic theory in t he solar wind are dis- 
cussed in detail in lHowes et al.l (|2008al ) . Energy spectra from 
the first nonlinear gyrokinetic simulations of Alfvenic tur- 
bulence at the scale of ion Larmor radius are modeled with 
striking agreement b y the cascade model used in this Letter 
l|Howes et alj|2008ch . 

A key simplification of the physics in the gyrokinetic 
limit is that the normalized eigenfrequency of the linear, 
collis ionless gyrokinetic dispersion relation (|Howes et al.l 
2006) is independent of fey, giving ZJ(k±, Pi, Ti/T e ) = 
co{k±, feti , pi, Ti/T e )/k\\VA\ the same holds for the damping 
rate, 7 = -y/kuVA- This enables the development of a one- 
dimensional model of the turbulent magnetic field energy (in 



velocity units) 



bt 



SB 2 i(k_ L )/4imim i , where <5B ± (k±)/8ir 



is the energy density of the magnetic field fluctuations per- 
pendicular to the mean field integrated over fcy . 

The model assumes a fully ionized plasma of protons 
and electrons with Maxwellian equilibrium velocity distribu- 
tions of temperature T s (in energy units). Plasma thermal 
velocities are taken to be non-relativistic, vu = \j2T s jm a 
c. At the scales of dissipation (of order the ion Larmor ra- 
dius), the mean magnetic field can be modeled as a straight, 
uniform field of magnitude Bo. The turbulence is assumed 
to be sub- Alfvenic, driven isotropically at a wave number fco 
with an amplitude satisfying the critical balance condition, 
and balanced (equal Alfven wave energy fluxes in opposite 
directions along the magnetic field). The model follows the 
magnetic energy in the Alfvenic cascade as it transitions 
from Alfven waves to kinetic Alfven waves at k±Pi ~ 1 
jHowes et al.ll2008d ; lHowes!l2008l ; ISchekochihin et al Jl2009h ■ 
The numerical algorithms used to calculate the linear gyroki- 
ne tic and Vlasov-M axw ell dispersion relatio ns are described 
in lQuataertl (|l99Sl ) and lHowes et all |2006t ). 

The continuity equation for the evolution of the mag- 
netic energy versus perpendicular wavenumber is given by 

M_ 
dt 



ge(fcx) 
dk_L 



+ S(k 1 



2 7 &L 



(1) 



1 Slightly different predictions for inertial range spectral slopes 
k ± 5 ^ or k^^ 2 arc unlikely to yield significant differences in the 
turbulent heating ratios predicted by this model. 



where the energy cascade rate is e(k±) = C 1 ~ 3//2 fcx^&i, the 
energy injection rate is S (non-zero only at the driving scale 
k± — feo), and the linear kinetic damping rate is 7. Assuming 
critical balance at all scales, the steady state solution for the 
energy cascade rate is given by 



j f kx n „3/2„ 7(fex) dk' x 
e(k ± )=e ex P {- / 2(7/ C a =^^- 



(2) 



where Ci and C2 are order unity dimensionless Kolmogorov 
constants and eo is the rate of energy input at ko- The solu- 
tion of e(fcj_) for parameters (Pi,Ti/T e ) is then used, along 
wit h the linear gyrokin etic damping rates due to each species 
7 S |Howes et al. l2006h . to calculate the spectrum of heating 
by species Q s (k±) = 2Cf //2 C r 2(7 s /cj)e(fcx)/fcx. After inte- 
grating over k±, we obtain the total ion-to-electron heat- 
ing rate due to the kinetic dissipation of the turbulent cas- 
cade, Qi/Q e (Pi, Tj/Te), the primary scientific result of this 
Letter. Note that iQuataert fc Gruzinovl (|l999l ) presented a 
similar model to estimate the fraction of electron heating for 
Ti/T e = 100 and 1 < ft < 10 3 , but did not account for the 
transition to a kinetic Alfven wave cascade at the scale of 
the ion Larmor radius k±pi ~ 1. 



2.1 Limitations of the Model 

First, the cascade at weakly collisional scales with k±pi < 1 
is assumed to consist entirely of Alfvenic fluctuations with 
Alfven Mach number less than unity. The energy in com- 
pressive modes, such as the collisionless manifestation of 
the MHD fast and slow magnetosonic waves, is neglected 
based on studies that demonstrate strong damping of these 
waves unde r weakly collisional conditions for som e parame- 
ter regimes l|Barnesll 19661 ; iFoote fc Kulsrudlll979T ); the appli- 
cability of this result over the entire range of parameters con- 
sidered here merits further study, beyond the scope of this 
work. Here we focus solely on the cascade of Alfven waves, 
which are incompressible to lowest order and remain un- 
damped down to the scale of the ion Larmor radius k± pi ~ 1 
(ISchekochihin et al . 2009). If significant turbulent energy ex- 
ists in the compressive wave modes, then its dissipation and 
the re sulting plasma h eating must be accounted for sepa- 
rately; [Quataert] l|l998h discusses this matter in more detail. 

Second, the entropy increase required to achieve irre- 
versible thermody namic heating can only be accomplished 
through collisions l|Howes et al.l 20061 ) . Wave-particle inter- 
actions damp the electromagnetic fluctuations comprising 
the turbulence and deposit the associated energy in to non- 
thermal features of the particle distribution function (|Howesl 
|200S| ; ISchekochihin et al] 120091) . It is assumed that the ac- 
tion of collisions on the non-thermal distributions of a given 
species lead to heating of that plasma species. 

Third, we do not explicitly account for the kinetic en- 
ergy of the Alfvenic fluctuations at k±pi < 1. The restoring 
force of Alfven waves is magnetic tension, so the energy in 
an Alfven wave is transferred back and forth between mag- 
netic energy and kinetic energy at the wave frequency; for 
an ensemble of many Alfven waves, the turbulent fluctua- 
tion energy is shared equally between magnetic and kinetic 
forms. As the turbulence transitions to scales k±pi > 1, the 
ions decouple from the turbulence, causing the kinetic en- 
ergy to become subdominant. It is assumed that the rapid 
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Figure 1. Contour plot ollog(Qi/Q e ) over the plane (ft, T;/T e ). 



transfer between kinetic and magnetic energy at the Alfven 
wave frequency enables both the kinetic and magnetic en- 
ergy at scales kj_pi < 1 to be channeled into the dominantly 
magnetic fluctuations in the kinetic Alfven wave regime 
k±pi > 1. Nonlinear simulations of kinetic turbulence will 
provide a crucial tool to test this assumption. 



3 RESULTS 

The key scientific result of this Letter, the calculation of 
ion-to-electron heating due to the kinetic dissipation of the 
Alfvenic turbulent cascade, Qi/Q s (/3i, Ti/T e ), is presented in 
Fig.[T]as a logarithmic contour plot over the parameter range 
0.01 j3i sC 100 and 0.2 s; Ti/T e sC 100. The values of the 
Kolmogorov constants used in this calculation, C\ = 1.96 
and C2 = 1.09, are derived from a fit of the cascade model 
predictions to nonlinear simulations of the tr ansition from 
the A lfven to the kinetic Alfven wave cascade dHowes et alj 
l2008d l. 

A salient feature of this result is that the heat- 
ing rate Qi/Qe is primarily a mon otonic function of Pi 
and is nearly independent of Ti/T e . IQuataert! (|l998l ) and 
IQuataert fc Gruzinovl l| 19991 ) provide an in-depth discussion 
of various kinetic damping mechanisms and their depen- 
dence on the plasma parameters; here we present only the 
key concepts necessary to explain this behavior. The ion 
damping peaks at k±pi ~ 1, while the electron damping 
peaks at k±pt ;$> 1 unless Ti/T e < m^/nrii, a condition un- 
likely to be relevant to astrophysical plasmas. Therefore, any 
energy that passes through the peak of the ion damping at 
k±pi ~ 1 will lead to electron heating. The damping onto 
ions via the Landau resonance is dominated by transit-time 
damping involving the interaction of the ion's magnetic mo- 
ment with the longitudinal magnetic field perturbation SB^ 
(the magnetic analog of Landau damping by the parallel 
electric field) l|Barneslll966l : IQuataert 1998j). As the mag- 
netic moment of the ions increases with increasing ft, so 
does the strength of the transit-time damping, leading to 
its dominant role in controlling the ion-to-electron heating 
Qi/Qe- The value of Ti/Te controls the scale at which the 
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Figure 2. Contour plot of the log(fcomaa;Pi) over the plane 
(0i,Ti/T e ), a necessary condition for the validity of the heating 
results presented here. 



electron damping occurs, but any energy passing through 
the range of ion damping will ultimately heat the electrons, 
leading to little effect on the ratio of total damping Qi/Qe- 

For values of T/Te < 0.2, the slow magnetosonic mode, 
often called the ion acoustic wave in this limit, is no longer 
heavily damped. In this case, the energy in this compres- 
sive mode may no longer be neglected. Further complicat- 
ing matters is the fact that the slow wave and Alfven wave 
may nonlin early exchange energy at t he scale of the ion Lar- 
mor radius (|Schekochihin et al.ll2009l) , so the ratio of ion-to- 
electron heating may no longer be accurately estimated us- 
ing the present model. Nonlinear kinetic simulations of the 
turbulence at the scale of the ion Larmor radius will yield 
valuable guidance to determine the turbulent heating ratio 
Qi/Q e in this regime. 

For the gyrokinetic treatment of Qi/Qe presented here 
to be applicable, the turbulence must be dissipated ki- 
netically (via the Landau resonance) before the frequen- 
cies rise to the point that the cyclotron resonance alters 
the dynamics. Since ui oc ku in the gyrokinetic limit, this 
constraint may be cast as a limitati on on the maximum 
isotropic driving wavenumber komax l|Howes et alj l2008al ) : 
for turbulence driven at wavenumbers below komax, fre- 
quencies remain sufficiently low that the cyclotron reso- 
nance never contributes and the gyrokinetic treatment is 
valid. The value of to m „ is determined using the following 
procedure: (1) using the critical balance condition to give 
fc|| = fc,!/ 3 k 2 / 3 Lu~ 1 ^ s (e/ eo) 1//3 , the maximum value of fcy/fco 
is determined from the steady-state solution e(fc±); (2) from 
this point (with fixed fcx, Pi, and Ti/T e ), the eigenfrequency 
using the linear Vlasov-Maxwell dispersion relation is solved 
as fey is increased until cyclotron resonant effects lead to 

[{iUvm, ~ U s k) 2 1 (w„m) 2 + (7„m ~ 7 9 fc) 2 / (7i>m) 2 ] 1/2 > 0.5; (3) 

the value of the parallel wavenumber k\\ max at this point 
is used to compute the corresponding ko max . A logarith- 
mic contour plot of the value of the normalized maximum 
isotropic driving wavenumber ko ma xPi over the (f3i,Ti/T e ) 
plane, presented in Fig.0 enables verification of the validity 
of the heating results for a turbulent astrophysical plasma 
with given values of /3i, Ti/T e , and fco- 
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Figure 3. Calculations of log(Qi/Q e ) vs. ft from the cascade 
model (dashe d), the fitted hea ting prescription eq. J3} (dotted), 
and cq. (5) of lQuataert] l|l99Sl ) (solid). 



3.1 Prescription for Turbulent Heating 

The second aim of this paper is to provide a prescription for 
the turbulent heating Qi/Qe presented in Fig.[l]that may be 
used as a sub-grid-scale model for heating in simulations of 
astrophysical plasma turbulence. The heating prescription, 
found by fitting the results for Qi/Qe in Fig. [I] is given by 



Qi/Qe = ci 2 

c§ + /3. 



mm 



m e T e 



(•3) 



where ci = 0.92, c 2 = 1.6/(Ti/T e ), c 3 = 18 + 5 log (Ti/T e ), 
and p = 2-0.2 log(T/T e ). A slightly better fit for Ti/Te < 1 
occurs using the coefficients ci = 1.2/(T/T e ) an d C3 = 18. 
It is in teresting to compare this result to eq. (5) of lQuataert! 
<|199SJ ). Qt/Qe = [(m I T I )/(m e T e )] 1/2 e- 1/ ' 3 ', which gives the 
predicted value of ion-to-electron heating by transit time 
damping in the limit k±pi < 1. In Fig. [3] log(Qi/Q e ) is 
plotted for Ti/T e — 0.2, 1, 100 for the three methods: the cas- 
cade model calculations (dash ed); the h e ating prescription, 
eq. © (dotted); and eq. (5) of lQuataert! dl99gf) (solid). One 
can see that, although the estimation by Quataertl (|l998l ) 
generally reproduces the behavior of Qi/Q e , significant de- 
viations at /3j ~ 1 are predicted by the cascade model calcu- 
lations, and the numerical results are much more accurately 
described by the heating prescription given by eq. ©. 



4 DISCUSSION 

We consider here the application of eq. © to studies of 
astrophysical accretion disks and to the interstellar medium. 

Recent two-fluid studies of the magnetorotational in- 
stability ( MRI) in accreti o n disk s using shearing box sim- 
ulations by ISharma et all (2007) have shown that the de- 
velopment of pressure anisotropics, and their saturation via 
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Figure 4. For a j3 = 10 plasma, the ion-to-electron heating ratio 
Qi/Qe predicted by eq J3]l (solid) and viscous heating Qv.il Qv,e 
arising from aniso tropic pressure due t o the background shear 
flow (dashed) from lSharma et all d2007T) ■ 



kinetic instabilities, leads to substantial viscous heating, rep- 
resenting a direct conversion of kinetic energy due to differ- 
ential rotation into heat at large scales. About 50% of the 
heating in the simulations occurred via this mechanism, with 
the ratio of ion-to-electron viscous heating given approxi- 
mately by Qv,i/Qv,e ~ 3(T;/T e ) 1//2 . Numerical grid-scale 
dissipation accounted for the remaining 50% of the turbu- 
lent energy loss; this loss was not added to the plasma heat- 
ing because it was unclear which species should receive this 
energy. By applying the predicted viscous electron heating 
to Sgr A* using a ID model for the radiative efficiency, the 
observed luminosity of L ~ 10 36 ergs s _1 implied an accre- 
tion rate of M ~ 10~ 7 to 10~ 8 M a yr" 1 , well be low the 
estimated Bondi accretion rate l|Sharma et al.ll2007l ). 

The 50% of the energy lost numerically in their simu- 
lations represents the energy transferred to small scales via 
turbulence, eventually to be damped by collisionless dissi- 
pation mechanisms. Here we estimate the partitioning of 
this heating by species using the turbulent heating model 
in eq. (3). The MRI drives turbulence on the scale height 
of the disk (|QuataertHl998T l. so we take the isotropic driv- 
ing wavenumber to be io ~ 1/H ~ 10 -8 km for Sgr A*. 
For an ion temperature of T, ~ 10 12 K, and a magnetic 



field strength of B ~ 30 G, the ion Larmor radi us is ap- 
pro ximately Pi ~ .4 km ijSchekochihin et al.ll2009l ). Follow- 
ing [Sharma|et^al] l|2007l ). we take the plasma Pi ~ 10 and 
Ti/T e ~ 10, so for a value of kop t ~ 10 -7 , we see from Fig. [5] 
that we are well below the maximum value komaxPi ~ 1 
required for the validity of eq. ©. 

In Fig. [4] we compare the turbulent heating ratio Qi/Q c 
from eq. Q for /3 = ftfl + Te/T) = 10 to the viscous heat- 
ing Qv,i/Qv,e from lSharma et all (|2007l ). The quite surpris- 
ing result is that, for a /3 = 10 plasma, the two distinct 
mechanisms yield a heating ratio that agrees within a factor 
of two over the entire range of T%/T e . Therefore, adding in 
the plasma heating due to the damping of the turbulence 
via kinetic mechanis ms does not change the conclusions of 
ISharma et all (|2007h regarding Sgr A* discussed above. 

Note that, in addition to the turbulent heating pre- 
scription presented here, a complete sub-grid model for tur- 
bulence would also require prescriptions for the sub-grid 
Reynolds and Maxwell stresses. The relative insensitivity 
of the measur ed stresses to the num erical resolution in the 
simulations of Sh arma et all (|2007l ). however, suggest that 
the unresolved fluctuations in those simulations do not con- 
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tribute substantially to the stresses. In this case, the heating 
prescription alone may be sufficient to determine the fate of 
the turbulent energy lost numerically. 

Turbulence in the interstellar medium can be probed by 
interpreting the electron density fluctuations as a passive 
tracer mixed by the ambient turbulence. Interstellar scin- 
tillation measurements of these density fluctuations in the 
diffuse ISM demonstrate an ou ter scale turbulent wave num- 
ber fco ~ 1/Lq ~ 10" 15 km ijArmstrong et alJll995l ). Al- 
though the ISM is an inhomogeneous plasma characterized 
by several different phases, it is the warm ionized compo- 
nent of the ISM that responds to the turbulent electro- 
magnetic fluctuations. Taking typical plasma parameters of 
n i = n e ~ 0.5 cm" 3 Tj ~ T e ~ 8000 K, and B ~ 2 x 1 0" 6 G 
[jSchekochihin et alJl2009l ). we find p t ~ 500 km and # - 3. 
For pi = 3 and Ti/T e — 1, the driving wavenumber of 
kopi ~ 10~ 12 easily satisfies the constraint in Fig. [2] Ap- 
plying eq. © , we predict an ion-to-electron heating ratio of 
Qi/Qe — 1 due to the dissipation of turbulence in the warm 
ionized component of the ISM. As seen in Fig. [3j the pre- 
dicted heating ratio is a relatively strong function of Pi, so 
for any particular volume of plasma, the heating may vary 
as Pi varies, but it is clear from the plot that in the range 
around Pi ~ 3 neither the ions nor the electrons are domi- 
nantly heated. Although there is a substantial cold, weakly 
ionized component of the ISM, collisional coupling with the 
warm ions is extremely weak at the dissipative scales, so it 
is unlikely to affect this prediction significantly. 



5 CONCLUSION 

This Letter presents a calculation of the ratio of ion to elec- 
tron heating Qi / Q e due to the kinetic dissipation of Alfvenic 
turbulence in ast rophysical plasmas b ased on the turbulent 
cascade model of Howes et all (|2008al ). The results demon- 
strate that the heating ratio Qi / Q e is primarily a monotonic 
function of Pi and has weak dependence on Ti/T e , as seen 
in Fig. [T] A fit of the cascade model calculations is given by 
eq. {3J, a form suitable for use as a sub-grid-scale model for 
turbulent heating in simulations of astrophysical plasmas, 
potentially enabling the determination of observational sig- 
natures caused by radiation emitted from the turbulently 
heated plasma. Future work aims to employ nonlinear ki- 
netic simulations of the dissipation of astrophysical turbu- 
lence to test these predictions of the turbulent plasma heat- 
ing. 
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